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Abstract. Recently, we have applied for the first time the angular momentum and 
isospin projected nuclear density functional theory to calculate the isospin-symmetry 
breaking (ISB) corrections to the superallowed /?-decay. With the calculated set of the 
ISB corrections we found | Kid| — 0.97447(23) for the leading element of the Cabibbo- 
Kobayashi-Maskawa matrix. This is in nice agreement with both the recent result of 
Towner and Hardy [Phys. Rev. C77, 025501 (2008)] and the central value deduced 
from the neutron decay. In this work we extend our calculations of the ISB corrections 
covering all superallowed transitions A, /'^ 0+, T 1, ^ A, = 0+, T = 1, + 1 
with Tz = —1,0 and A ranging from 10 to 74. 
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1. Introduction 

The isospin symmetry [H [2] in a nuclear medium is only weakly broken reflecting 
the relative weakness of isospin-breaking as compared to the isospin-conserving part 
of the nucleon-nucleon interaction. Hence, the related isotopic spin quantum number, 
T, albeit approximate, remains very useful for labeling nuclear states and understanding 
selection rules on different types of nuclear reactions (see Ref. [3] and refs. therein). In 
particular, for the Fermi (vector) and Gamow- Teller (axial) /3-decay one has AT = 
and AT = 0,±1, respectively, with the exception of T^ = O"*" — ?■ T^ = 0"'' transitions 
which are strictly forbidden for the Gamow- Teller process [4j. 

Among pure Fermi transitions of particular importance are superallowed transitions 
between the isobaric analogue states I'^ = 0+,T = 1,T2 — )■ T" = 0~^,T = 1,T2 + 1 in 
N Z nuclei. These data are used for precision tests of the conserved vector current 
(CVC) hypothesis i.e. the independence of the vector current on the nuclear medium. 
With the CVC hypothesis being verified, they serve as the most precise source of the V^d 
element of the Cabibbo-Kobayashi-Maskawa (CKM) matrix which is a key ingredient 
in investigating its unitarity. 

The CVC hypothesis is verified by investigating nucleus-independence of the J^t- 
values: 

n = ft{l + 5'^)il + (5ns - Sc) = 2Gl(l + Al) ^ 

where K/{?icf = 2Ti^n\Yi2l {m^c^f = 8120.2787(11) x 10"^° GeV'^s is a universal 
constant and Gy stands for the vector coupling constant for the semileptonic 
weak interaction. The J-'t-values include empirical reduced life-times ft corrected, 
theoretically, for radiative processes and isospin-symmetry breaking. The radiative 
corrections are routinely divided into the following: a nucleus-independent part A^ = 
2.361(38)% [5], a transition-dependent (Z-dependent) but nuclear-structure-independent 
part [HI E] and a nuclear-structure-dependent part (5ns [3 E] . The isospin-symmetry- 
breaking correction, Sq, is defined through the following nuclear matrix element: 

|(/ = 0,T ^ 1,T, = ±1|T±|J = 0,T ^ 1,T, = 0)|2 = 2(1 - 6c), (2) 

where T± are the raising and lowering bare isospin operators, respectively. 

Application of the superallowed /3-decay in testing the three-generations-quark 
Standard Model of elementary particles requires both high-accuracy empirical /t-values 
and a high-quality theory. The aim of this work is to present new results on the 6c 
corrections obtained using angular momentum and isospin projected density functional 
theory (DFT) recently developed by our group P El HOl [IH [121 |13l E] . The model will 
be introduced in Sec. [2l the results on 6c will be presented in Sec. [3] and the paper will 
be briefly summarized in Sec. HI 
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2. The model 

The degree of isospin mixing in atomic nuclei (isospin impurity) is, predominantly, a 
result of subtle interplay between the short-range isospin-symmetry-conserving strong 
interaction and the long-range isospin-symmetry-breaking Coulomb force. The Coulomb 
force polarizes wave functions of all participating protons and in turn, neutrons, creating 
a number of conceptual and technical difficulties within the nuclear shell-model (SM) 
having rather profound consequences for calculations of 6c corrections [151 IIBl [13 El HE] . 
The SM can be used to compute only a part of the correction related to configuration 
mixing 6ci- The second part, 6c2, related to the radial mismatch of the single-particle 
wave functions, must be calculated independently using a mean-field. Both corrections 
are treated as additive: 6c = ^ci + ^C2- The consequences of this rather artificial division 
include, for example, a necessity of using effective isospin operators which violate the 
SU(2) commutation rules. This problem was recently noticed and discussed extensively 
by Miller and Schwenk [19] who, however, didn't provide any quantitative estimate of 
the impact of these deficiencies on the results by Towner and Hardy ^ (TH) whose 
calculations set the standard in this field. 

In contrast, Hartree-Fock (HF) and DFT are free from these specific problems. 
Here all nucleons participate on equal footing and a balance between long- and short- 
range effects is treated in a self-consistent manner. It is well known, however, that these 
approaches break spontaneously fundamental nuclear symmetries including rotational 
and isospin symmetry. In fact, the isospin symmetry is violated both explicitly, by 
virtue of charge-dependent interactions, and spontaneously which leads to unphysical 
quenching of true isospin mixing. 

In order to avoid spontaneous mixing and compute matrix element ([2]) in a fully 
quantum mechanical way using bare isospin operators we have recently developed the 
angular-momentum and isospin projected scheme on top of the Skyrme-DFT approach. 
The approach is based on rediagonalization of the entire Hamiltonian including the 
Coulomb interaction in good angular momentum and good isospin basis 



the isospin and angular momentum projection operators (see [20] for further details). 

For the Fermi matrix element the normalized bra and ket states are calculated 
in the following way. The |/ = 0, T f» 1, = ±1) state in the even-even nucleus is 
projected of the self-consistent Slater determinant, {ip), representing the ground state 
in this nucleus: 



The state \ip) is unambiguously defined. The \I = 0,T ^ l,Tj = 0), on the other 
hand, is projected of the self-consistent Slater determinant, \ip), representing the so- 
called antialigned configuration \{p) = |P (g) tt) (or |z/ (g) tt)), selected by placing the odd 





(4) 
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neutron and odd proton in the lowest available time-reversed (or signature-reversed) 
single-particle orbits: 

|/ = 0,T ^ 1,T, = 0) = E c^^^P^oPoflv)- (5) 

T>0 

The selected single-particle configuration \u^it) manifestly breaks the isospin symmetry 
(see Fig. 1 in Ref. [H]). This is essentially the only way to reach the \T 1,1 = 0) 
states in odd-odd N = Z nuclei. Indeed, only the T = states in N = Z nuclei can be 
represented by a single Slater determinant. The final expression for the nuclear matrix 
element (E]) will be given in our forthcoming paper. 

The two major drawbacks of the model in its present formulation include (i) the 
lack of pairing-correlations and (ii) the use of an old fashioned and low-quality SV 
parametrization of the Skyrme force. The latter deficiency pertains to the angular 
momentum projection which is known to be ill-defined for density-dependent modern 
Skyrme and Gogny energy density functional (EDF) [211 1221 [231 IMl [25]. At present, 
the SV interaction augmented by a tensor term is essentially the only available Skyrme 
interaction originating from the true Hamiltonian which can be safely used in connection 
with the angular momentum projection without any further regularization [25j. It is 
worth mentioning here that the isospin impurities calculated using modern Skyrme forces 
in the isosp in-projected variant of our model [101 [12] are consistent with the recent data 
extracted from the Giant Dipole Resonance decay studies in ^°Zr [26] and the isospin- 
forbidden El decay in ^^Ge [27]. This indicates that our model is in principle capable of 
quantitatively capturing the amount of isospin mixing that is important in the context 
of making reliable calculations of 6c in spite of the fact that 5c is mostly sensitive to 
the difference between the isospin impurities of parent and daughter nuclei [281 129]. 

3. The results 

The isospin-breaking corrections Sq were computed by different groups and with diverse 
nuclear models [30l [6l [3ll [321 129]. Our calculations were published recently in Ref. [H] 
and compared to the TH [6]. Comparison of 6c s shows that although individual 
values differ, both sets of calculations follow a similar trend with increasing A. The 
differences between individual values are stronger in light nuclei and are traced back to 
poor spectroscopic quality of the SV parametrization. It is interesting to note that our 
results are considerably larger as compared to the results of Ref. [32] which are based 
on the relativistic Hartree (RH) plus random phase approximation (RPA) formalism. 

The average value of the nucleus-independent reduced hfe-time ([H) calculated for 
twelve out of thirteen high-precision superallowed /3-decays (excluding the ^^K— i-'^^Ar 
transition) equals J-')f:=3070.4(9)s, which is consistent with the CVC hypothesis. This 
value was obtained using Gaussian-distribution-weighted formula to conform with the 
standards of the Particle Data Group (PDG) and using our adopted 6c listed in Table 1 
of Ref. [11]. In the calculations we used the radiative corrections and /t- values taken 
from Ref. [6] and [18], respectively. The calculated J-'t leads to |Vud| = 0.97447(23) 
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Figure 1. Differences between theoretical and experimental Q-values for the 
superallowed ;3-decays r = 0+,T = T = 0+ ,T = 1,T^ + I versus the 

atomic number A. Black and white symbols represent transitions in mirror-symmetric 
pairs corresponding to = — 1 and = 0, respectively. 



which is in nice agreement with the TH result \V^^ \ = 0.97418(26) and a central value 
deduced from the neutron decay |V^ud''l = 0.9746(19). Combining our |Vud| with the 
results of = 0.2252(9) and |Kb| = 0.00389(44) adopted from the most recent PDG 
compilation [33] leads to iKidP + iKisP + |Kib| = 1.00031(61). This result implies that 
the CKM unitarity inferred from its first row is fulfilled with a 0.1% precision. 

The confidence level (CL) of our calculations can be assessed be performing the test 
proposed recently in Ref. [H]. The test indicates that the CL of our model is lower than 
the CL of other models analyzed in Ref. [I8] . It should be stressed, however, that the low 
CL of our model results primarily from a very small 6c for a single transition in A=62. 
This seems to be a part of deeper problem faced not only by our model but, in fact, 
all models involving a self-consistent mean-field. Indeed, these models have problems in 
reproducing quantitatively the rapid increase in Sc expected to occur already at A ~ 62, 
see [13. In our model the increase in Sc takes place at A = 70, similar to the RH plus 
RPA calculations of Ref. [32] ■ To lesser extent this problem is also seen in the SM plus 
HF calculations of Ref. [HI [3l]. 

The low confidence level of our results reflects poor spectroscopic quality (ordering 
and energy of the single-particle levels) of the SV EDF. Global characteristics of the 
functional are also unsatisfactory. In particular, it reproduces the nuclear binding 
energies of even-even = ±1 partners (after projection) with accuracy of order of 
1% — 2% relative to the data. As already mentioned 6c probes differences between the 
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Table 1. Calculated values of Sq for the superallowed /3-decays — 0+,r = IjT^ — > 
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partners and these quantities i.e. the relative binding energies between the isobaric 
analogue states in the /3-decay partners, the Q-values, are captured quite well (see 
Fig. [1]). Nevertheless, even the Q- values show, in some cases, relatively large deviations 
from the data and rather unacceptably large fluctuations versus A, in particular, in light 
nuclei. Part of these deviations reflects the well recognized deflciency of conventional 
Skyrme-DFT models known in the literature under the name of the Nolen-Schiffer 
anomaly, see [351 ESI EZ] . In the calculations it manifests itself as a discontinuity in the 
calculated relative Q- values between the = — 1 — )■ = and = — )■ = 1 
transitions as shown in Fig. [H Indeed, the flgure shows that theoretical Q-values 
tend to underestimate the experimental data and that the effect is larger in the 
Tj; = — 1 — )■ = cases as compared to the = — )■ Tj = 1 transitions. At 
present it is neither clear how to cure this anomaly nor how to estimate its influence on 
the isospin mixing. 

In order to check for the stability of our predictions we have carried out systematic 
calculations of 5c foi' all = —1 — t- = and = — ?■ T^, = 1 transitions 
from A = 10 to A = 74 being aware that most of these cases will probably never be 
measured. In the calculations we used = 10 and = 12 spherical harmonic oscillator 
shells for A < 34 and A > 38, respectively. Moreover, the anti-aligned states in N = Z 
odd-odd nuclei were calculated assuming signature-symmetry. The present results were 
calculated under different conditions and are therefore different than our preferred values 
given in Ref. [13] • They are collected in Tab. [Hand visualized in Fig. O 

The calculations conflrm earlier results concerning large mixing in A = 18 and 
pathologically large mixing in A = 38. Unusually large mixing is also obtained in 
A = 58. Note that these difficulties are seen in both the = — 1 — t- = and the 
T2 = — 7- = 1 cases. In general, the values of 6c in mirror-symmetric transitions 
Tz = —1 ^ Tz = and = — )■ = 1 follow a similar trend with increasing A 
but the individual values are different. In some cases the differences are sizeable. It 
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Figure 2. Isospin-mixing corrections Sq to the superallowed /3-decays I"^ = 0+,T = 
l,Tz — >/ = 0+,T = l,Tz + 1 versus the atomic number A. Black and white symbols 
represent transitions in mirror-symmetric pairs corresponding to Tz = — 1 and Tz — 0, 
respectively. 

is interesting to observe that the largest differences between Sq in mirror-symmetric 
pairs occur in A = 62 and A = Q6 i.e. exactly in the region where the transition from 
small to large values of 6c is predicted to occur in the SM plus Woods-Saxon potential 
calculations |6l HE]. This result indicate that configuration mixing around A = 62 is 
very fragile in the self-consistent calculations. The effect may be sensitive to various 
characteristics of the underlying EDF and/or to missing correlations and requires further 
study. 

4. Summary and conclusions 

We have extended the isospin and angular momentum projected DFT calculations of 
the isospin breaking corrections 6c to all theoretically possible superallowed /3-decays in 
nuclei ranging from A = 10 to A = 74 in order to assess the stability of our predictions. 
The calculations reveal that values of 6c in mirror-symmetric decays Tj = — 1 — ?• = 
and Tz = ^ Tz = 1 follow a similar trend versus the atomic number A. The individual 
corrections in mirror-symmetric decays are shown to be different, in some cases sizeably, 
reflecting differences in configuration mixing in = ±1 nuclei. The largest differences 
were found for cases A = 62, 66 i.e. in the region where the SM plus Woods-Saxon 
potential calculations [HI [18] predict rapid transition from small to large values of ^c- 
This transition is delayed in the self-consistent calculations. 
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